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Abstract. Taking the XXZ chain as the main example, we give a review of an 
algebraic representation of correlation functions in integrable spin chains obtained 

• recently. We rewrite the previous formulas in a form which works equally well 

f^*) , for the physically interesting homogeneous chains. We discuss also the case of 

■ quantum group invariant operators and generalization to the XYZ chain. 
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£v] ■ 1. Introduction 



fSJ ! The investigation of integrable spin chains has a long history since Bethe's work 
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PQ, in which the Bethe Ansatz method was invented. It was only a start, and later 
was followed by a line of new ideas and concepts such as commuting transfer ma- 
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trices, the Yang-Baxter equation, the quantum inverse scattering method, quantum 
groups and the quantum KZ equation. In a series of papers 121121 SIS], we studied an 
algebraic formula for the correlation functions in the infinite XXX, XXZ and XYZ 
spin chains. Our method is a synthesis of those mentioned above. 

The study of correlation functions has been a highlight in the researches of these 
spin chains. In the early days the only knowledge was the nearest neighbor cor- 
relator, written in terms of log 2 for the XXX model. It was a big surprise when 
Takahashi f2H] found ((3) in the next-nearest correlator, where £(s) is the Riemann 
zeta function. In ^T] ^1 EI], the quantum vertex operators in the representation 
theory of the quantum affine slz algebra were used to obtain multiple integral for- 
mulas for the general correlation functions of the XXZ model. Kitanine, Maillet, 
Slavnov and Terras rederived and further generalized these integral formulas to in- 
clude magnetic field and time f2J El (see PB] for a review). Study of the finite 
temperature case has also been launched recently by Gohmann, Klumper and Seel 
|lUj . and progress has been made in the calculation of long distance asymptotics of 
some correlators by the Lyon group jTTj and Korepin, Lukyanov, Nishiyama and 
Shiroishi [T^j. However it was not immediately understood why £(3) appears in the 
next-nearest correlators. 

In (HI U\ Boos and Korepin explicitly performed the multiple integrals for the next- 
nearest case and beyond, and again found odd integer values of C( s )- Further exact 
results including the XXZ chain have been obtained by Kato, Nishiyama, Sakai, 
Sato, Shiroishi and Takahashi |2U 1213 I2IH !2Zj • I n |Sj Boos, Korepin and Smirnov 
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studied the inhomogeneous correlation functions for the XXX model, and arrived 
at a conjecture on the algebraic structure for the general correlation functions: in 
brief, one transcendental function is enough to describe all of them. In the limit of 
the homogeneous chain, the Taylor series expansion of this function produces the 
special values of £(s) as well as log 2. 

In 013111, we proved the conjecture by giving an algebraic formula, and obtained 
similar results in the XXZ and XYZ models. The number of transcendental functions 
increases to two and three, respectively, as the number of parameters in the models 
increases. The main idea in the proof was the use of the reduced quantum KZ 
equation, and the main ingredient in the algebraic formula was the transfer matrix 
defined via an auxiliary space of non-integer dimensions. 

The algebraic formulas presented in these papers had some deficiencies: the beauty 
of the formula was marred by a chip on the edge of a comb. The relevant transfer 
matrices are 'incomplete', in that they act on the tensor product where two spaces 
are omitted. Also the formula for the inhomogeneous model consists of a sum of 
terms, which have poles when one tries to take the homogeneous limit. They cancel 
each other only after the summation. In [3] , these spots were cleaned up in the XXX 
model. 

In the present paper, we give the algebraic formula for the density matrix in a 
transparent form, not only in the XXX model but also in the XXZ and the XYZ 
models. We use the infinite XXZ chain as the main object: 



The density matrix p n belongs to the space End((C 2 ) (gln )* dual to the space of local 
operators End((C 2 ) lg)n ). We consider the space (C 2 )® n as the subchain of the entire 
infinite chain on which the XXZ Hamiltonian acts. It has the defining property 



Here the right hand side is the ground state average of the operator e End((C 2 ) (gin ). 
We will give the formula for p n in the form 



x (u;i(/ii >2 )X[ |n (/zi, /z 2 ) + w 2 (Aii,2)3^, n (^i, 

Here, k is a constant, f?(/ii i2 ) is a 4 x 4 matrix depending on /z 12 = pi — /^ 2 , an d 
Ui (i = 1,2) are certain transcendental functions. The operator T*(/i) is given in 
terms of the L operator L(p) G £/ ? (st 2 ) <8> End(C 2 ) as the monodromy matrix in the 
adjoint action: 



(1.1) 




p n (0) = (vac|0|vac). 





■ ■ ■ LM^OLM ■ ■ ■ Li(fi) £ U q {sh) <g> End((C 2 )^). 
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The deformation parameter q = e m " and the anisotropy parameter A in (jl.lj) are 



The operators X* n (/j,i, H2) {i = 1,2) is obtained from the monodromy matrix 



where Tr^ denotes the (i-dimensional trace. 

The above formula is in the homogeneous case, and the integrand has poles at 
Hi = (i = 1,2). The integral means taking residues at these poles. A similar 
formula is also given in the inhomogeneous case where the Hamiltonian is replaced 
with the transfer matrix for the inhomogeneous six vertex model with the spectral 
parameters Ai,...,A n associated with the tensor components of (C 2 )® n . For the 
details, see Theorem 13.11 and ()4.2j) - (j4.3|) . In this case, the integrand has poles at 
Hi = Xj (i = 1, 2; j = 1, . . . ,n). Taking residues at these poles we get the formula 
obtained in the previous paper 

For a general local operator 0, we need two functions u>i (i = 1, 2) to express its 
expected value. However, if is invariant under the action of [/^(sUi), the formula 
simplifies, and we need only u\. This case is related to the spin chain with an 
open boundary condition given by the Pasquier-Saleur Hamiltonian [22] • The XXZ 
Hamiltonian with periodic boundary condition corresponds to the CFT with the 
central charge c = 1. In contrast, the Pasquier-Saleur Hamiltonian corresponds to 
the CFT with c = 1 — 6z/ 2 /(l — v). The above property of the invariant operators 
was conjectured in j^j. We give a proof to this conjecture. 

We also give a formula similar to (jl.2j) . (J1.3|) for the XYZ model. 

The paper is organized as follows. In Section |21 the density matrix is defined. In 
Section El an algebraic formula of the operator Q n , which is dual to fi*, is given. In 
Section 01 the algebraic formula is written in an alternative form. In Section the 
formula for the invariant operators are given. In Section |U1 the formula for the XYZ 
model is given. 

The text is followed by three appendices. In Appendix [21 we give the derivation 
of the new formula for Q n . In Appendix El we make a comparison between different 
conventions used in this paper and in the book 12]. In Appendix O formulas for 
the normalization factors are gathered for the XXZ and the XYZ models. 





2. Density matrix for the XXZ chain 



Consider the XXZ Hamiltonian 



(2.1) 




where a a (a = 1, 2, 3) are the Pauli matrices and 



A = COS7TI/ 
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is a real parameter. We consider the two regimes, the massive regime A > 1, 
v G iM>o, and the massless regime |A|<1, 0<i/<1. 

Take a sub-interval of the lattice consisting of sites 1, . . . , n, where n is a positive 
integer. Let (E e ^)j denote the matrix unit (S ae 6be) a b= ± acting on the site j. By 
a density matrix, we mean the one whose entries are the ground state averages of 
products of the (E e ^)j's, 

(2.2) p n (A 1; ...,A n ) 

= Yl A 1 ,...,A 7I (vac|(£ ?1 , £1 )i • • • (E Zn}en ) n \vac)x 1 ,...,x n {E ei ,e 1 )i • • • (£e„,ejn- 

fl,-,«n 



Here we consider the model with inhomogeneities Ai, • • • , A n attached to each site. 
More precisely, we mean the following. 

Let V = C 2 be the two dimensional vector space with basis v + ,v_. Throughout 
this paper, we set 

(2.3) q = e™. 

Denote the standard trigonometric R matrix by 

P(A) 



(2.4) 



R(X) 



r(A) 



[A + l] r(A) ' 

/[A+l] 

[A] 

1 

\ 





1 

[A] 










\ 



e End(V® V). 



[X + 1]J 



Here the entries are arranged in the order (++), (H — ), ( — h), ( ), and 



[A] 



q — q 
q-q- 



The factor p(A) = p(A, 2) will be given later (see fH|) and ffTTTjl . (f(T2jl ). Introduce 
an auxiliary space V a — V with spectral parameter A, and denote by R a j the 
R matrix acting on V a ® Vj. Using (|2.4j) . we consider the transfer matrix of the 
inhomogeneous six vertex model 

(2.5) tx Va {R a>L {X - Ax) ■ • -i4 jn (A - An) ■ ■ ■ R a ,i^ ~ Ai) • ■ -Ra,-^ ~ 
which acts on the tensor product 

V-l ® ■ ■ ■ ® Vl ® ■ • • ® v n ® - • • ® V L . 

With each Vj we associate a spectral parameter Xj, assuming for defmiteness that 
Xj = for j < or j > n + 1. Let |vac)^... iAn denote the eigenvector of (|2.5jl corre- 
sponding to the lowest eigenvalue. We denote the dual eigenvector by A 1 ,...S(vac|, 

normalized so that x 1 ,...,^(vac|vac) A f'^. iAn = 1. The vacuum expectation value in 
()2.2|) is defined to be the thermodynamic limit 

A 1 ,...,A„(vac|(£ ?li£1 )i • ■ • (£ ?n , £ J„|vac) Al ,..., An 
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For an arbitrary local operator G End(V®") , we have 
( 2 -6) A lv ..,A„(vac|0|vac) Ali ... iAn = tr y ®n(Op n ). 

Our aim is to give an algebraic representation for the density matrix p n . 

3. Algebraic formula 

The density matrix p n is an operator on V® n . To present the result, let us pass 
from operators to vectors in V® 2n . We number the spaces as 

(3.1) Vi <2 • • • <2 V n <2 V* <2 • • • ® V T . 

We use the following convention for the indices: for example, if u = u ' ® u " ■> v = 
v' <2> v" are vectors in V <2> V, then we write 

M 1,T W 2,2 = «' <8> ® <8> U" e Vi <2 ® ^2 ® ^T- 

Similarly, we indicate by suffix the tensor components on which operators act non- 
trivially. 

Introduce a function h n with values in (j3.1j) : 

n 

(3.2) ^(Ax, ...,X n ) = J2 n(-e;)(vaci(£-e l!ei )i • • • (£-, ft , e J n |vac> 

3=1 

xu £1 (2 • • • <g> u £n (8) v tn <2 ■ ■ ■ <2 vz x . 

In Section we discuss more about the transition from p n to h n , and vice versa. 
We mention here only that spectral parameters Xj, Xj + 1 are attached to the spaces 
Vj and Vj, respectively. 

The function h n is known to satisfy the following system of equations: 

(3.3) h n (. . . , Afc + i, Afc, • • •) = Rk,k+i(^k,k+i)Rk^i ; k{^k+i,k)h n (. . . , A&, A^+i, . . .) , 

(3.4) h n (Xi — 1, A2, • • • , A n ) = A n (Xi, . . . , X n )h n (Xi, A2, • • • , A n ) , 

(3.5) T l J h n (Xi, . . . , A n ) = -s 1 j/i n _i(A2, • • • , A n )2,..., n ,n,...,2 > 

•Pn,nAn(Ai, • • • , A n ) = -S„^/l n _i(Ai, . . . , A n _i) 1) ... jn -l,^T,...,T ■ 

Here the notation is as follows. We set Ajj = Xi — Xj, R = P R with P being the 
transposition, 

*- = \(I-P) 
is the projection onto Cs where s denotes the vector 

(3.6) s = v + ® V- — V- ® v + E V ® V, 
and 

Ai(Al, • • • , A n ) = ( — l) n P 1 J-R li 2(Ai i 2 — 1) • ■ ■ -Rl,n(Al, n ~ l)-Rl,n(Al, n ) • • • #1,2^1,2) • 

^There is an erratum in [3]; the right hand side of the formula seven lines below (13.1) should 
read ]Xj =1 {~ej)(va,c\(E^ uei )i ■ ■ ■ (£L ?nilE J n |vac). 
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We call (|3.3J1 - (|3.5|) reduced qKZ (rqKZ) equations. We are going to construct 
a solution of these equations in a certain specific form. For that purpose, we will 
need three ingredients: monodromy matrix, trace functional Tr A , and transcendental 
functions u),u 2 (X, u). Let us explain them. 

Let E,F,H be the standard generators of U q (sl 2 ). We consider the L operator 
with auxiliary space of dimension d, 

(3.7) L W = 

where 

£(X) = ( [X t^f ] n F ^Z H ) 6 U q (sl 2 ) <g> End(V). 

\ q 2 E [X + ±f-\J 

The normalization factor p(X, d) is chosen to satisfy 

L(X)L(— A) = 1 (g) iy (unitarity relation), 

<x 2 L(A)V 2 = — L(— 1 — A) (crossing symmetry). 

Here L(Xf is the transposed matrix with respect to End(V). For the explicit formula 
of p(X,d), see (|C.1|) . (jC.2|) . We have, in particular, 

(3g) p{\d) p(X -l,d) 1 

{ ' } [A + f][A + f-l] [A-|][A + |]' 

We define the monodromy matrix T n (X) — T n (X\Xi, . . . , X n ) by 

(3.9) T n (X) = Li(A - Ai - 1) • • • L n (X -X n - l)L n (X - X n ) ■■■L 1 (X-X 1 ). 
The trace functional Tr A is the composition map 

(3.10) Tr A : U q (sl 2 ) - U q (sl 2 )/[U q (sl 2 ), U q {sl 2 )} - AC[C, C" 1 ] © C[C, C 1 }, 

where ( = q x . The first map is the canonical map, and the second is defined by 
setting for any m G Z 

mH ^ _ J [mA]/[m] if m + 0; 

if m = 0, 



Tr A (g* 

and for any x G t/^sfe) 




Tr A (Ca;) = 9 + 9 Tr A (x), 
(g-g- 1 ) 2 



where C is a central element given by 



C =-f t ^ + ER 



An equivalent way of defining Tr A x for x G U q (sl 2 ) is as follows. It is the unique 
element of AC[C,C _1 ] © C[C,C _1 ] suc h that, for all k + 1 dimensional irreducible 
representation n^ k ' : U q (sl 2 ) — > End(C fe+1 ) we have 

(Tr A x) |a=A;+i = tr y(fc) 7r (A:) (x) (A; G Z> ). 
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With this definition of Tr^, the 'trace' of the monodromy matrix has a unique 
decomposition 



(3.11) Tr W2 T n 



2 



— Xi tTl (fJ,i, /i 2 |Ai, . . . , A n ) — /il,2^2,n(A t l) A*2 1 ^1 s ■ ■ ■ j ^n)j 

where X ijn ( / u 1 , /i 2 ) = X i) „( / u 1 , /i 2 | Ai, • • • , A n ) (i = 1,2) are matrices whose entries 
are rational functions in the variables g Ml , g M2 , q Xl , ■ ■ ■ ,q Xn . Note that, with the 
substitution A = Ml + M2 — Xj, d = /i lj2 , the right hand side of (J3.8j) becomes 

( 3 - 12 ) -r A n- 

Finally, define the functions Ui(X, v) {% = 1, 2) 2 by 
(3.13) Kdlogv^A, i/) = ll>i(A, z/)rf(Az/) + c<j 2 (A, u)du, 

where 

/[A-l]\ 1/4 sin to; 

^(A, 1/) = p(A) [j— — T , k- 



.[A + iy ' To/ ' 

and v is given in ()2.3|) . The function p(A) = p(A, 2) depends also on z/ and is defined 
in ((C.IJ) . (|C.2|) in each regime. 

We are now in a position to state the algebraic formula. Set 



(3.14) Sn = H 

Theorem 3.1. The following formula gives a solution of the rqKZ equations (13. 3|) - 
iB: 

(3.15) h n (\ u --- ,X n ) = ^e n ^-^s n , 
where 

'-l) n f f dfii d\i 2 



(3.16) O n (Ai, • • • , A n ) . . 

2k 2 y y 27u 27n 

x (^i(A*i,2)X lin (/ii,// 2 ) + ^ 2 ( / ui i2 )X 2i „(/ii,/i 2 )) Tr 2j2 (T n (/ii) ®T n (^ 2 ) ■ B(^ 1)2 )) 

and 

/0 \ 

g/^ _|_ q-n _q _ g-i 



— q — q 1 q^ + q ^ 



(3.17) = - 7 M r 

1 ' W 2 |>* - 1] [>* H- 1] 

V °/ 

Hereui(\, v) are given in ()3.13|) . X ajn are defined by 1)3.11)1 . Tr 2j2 (x®|/) = Tr 2 (x)Tr 2 (?/), 
and t/ie matrix B(fi) acts on the auxiliary space C 2 <S>C 2 . TTie integral in fi{ (i — 1,2) 
means taking residues at the poles fj^ — Ai, • • • , A n which result from \3.1fy) . 
In the massive regime, h n coincides with the vector form of the density matrix. 

2 Our cl>i, cl>2 here are denoted u> and u> in )S]. 
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Conjecturally the same formula gives also the density matrix in the massless 
regime as well. 

In the above formula, we abbreviated all the variables u, X±, etc., other than H\,/j,2- 
This formula was given earlier in |3] in a different form. Since the integrand has no 
pole at Aj = Xj, the present formula is equally valid in the homogeneous chain where 
Ai = ■ • • = X n = 0. We show the equivalence of the two formulas in Appendix 1X1 

There is another representation for the operator Q n using the 'invariant' trace 



(3.18) 



Tt g x (A) = Ti x (q- H A). 



Define ,X| n and by 



(3.19) 



Tr q T 
- L Vi2- tn 



Mi + A*2 



X? n (//i, M2IA1, 



(3.20) 



A n ) — A i l,2X| n (/ii, ll2\X\ 

/0 \ 



An), 



[//-!][//+ 1] 



q 



q 



.qn q 



\ 



Then fl n can also be written as (|3.16|) . with X a>n (fjLi, ^2) and B(fi) replaced by 
^0,71(^1)^2) an d B q (fi), respectively. 

The existence of this second representation is a peculiar feature of the XXZ model 
which has analogs neither in the XXX nor in the XYZ models. 

Remark 3.2. The choice of the operator B(fi) in ()3.17|) is not unique. For example, 
there is a freedom of adding identity to B(/i) (see Lemma \A.5\) . 



4. An alternative representation 

In this section we return from vectors in V® 2n to operators on V® n . Let us recall 
some generalities concerning the action of quantum groups on these spaces. 
Denote by tt\ : [/'(sk) — > U q {sii) the evaluation homomorphism 

7r A (e )=g A F, 7r A (/o) = q~ X E, n x (^) = q^\ 
7T X (e 1 ) = q x E, n x (f 1 )=q' x F, 7r x (q ±h ^) = q± H ' 2 . 

For the representation ir^ of U q (sl2), we set n x = tt^ o 7t\. We use the coproduct 

A(e<) = ei ® 1 + q hl ® et, 
A(/i) = /i^g-^ + l®/* 

to define the action of [/'(sfe) on a tensor product of representations. 

Quite generally, for a finite dimensional [/'(sfe) module its dual vector space 
W* has two module structures defined via the antipode S as 

(xu, v) = (u, S ±1 (x)w) (x e C/,(sI 2 ), w eW*,veW). 
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Denote these structures by W* s±1 . We have canonical isomorphisms 

(w**)**' 1 ~ w, (Wi ® w 2 )** ~ wj* ® ^r* 

for = S ,±1 . The canonical pairing jy* 5 ® PU — * C is [/'(s^) -linear. We regard 
^(sb) -module via 

End(lU) ~ W ® TU* 5 

Using the trace tr^(AB), End(TU) may be identified with its dual space. The 
induced dual module structure becomes 

End(TU)* 5 " 1 ~ (W ® IU* S )* 5_1 ~ ® TU* S_1 . 

We are mainly concerned with the 2 dimensional module V where the generators 
E,F,H act in the basis t>+, i>_ as 

E = (0 0) ' F = (l 0) ' H = (0 -l) • 

Use the letter V(A) to indicate the evaluation module structure tt^ on V . We have 
then an isomorphism of U'Js I2 ) -modules 

K(A)* S±1 ~y(A=Fl), < 1- et;_ £ , 

where (v*,v e >) = 8 e>€ >. In particular, the identity operator I v G End(K) corresponds 
to s G U <g> V given in ([3.6)1 . 

We started from the tensor product 

S n = V 1 ®---®V n , V j = V(X J ), 

corresponding to the finite interval 1, . . . , n on the lattice. Our space of local oper- 
ators is the t/g(s[ 2 ) -module 

(4.1) L n = End(S n ) 

~ K(Ai) ® • • • ® U(A n ) ® U(A n - 1) ® • • • ® U(Ai - 1) 
on which x G [/'(s^) operates by the adjoint action 

ada;(0) = J2 x 'i 0S ( x ") (° G £ ™)> 
where A(x) = ® x i ■ I n contrast, density matrix belongs to the dual module 
L* = End(S n )* s_1 

~ V{X 1 ) ® • • • ® V(X n ) ® U(A n + 1) ® • • • ® V(Xi + 1). 

The action of x G ^/'(s^) is, in the same notation as above, 
ad'x(CT) =Y,<0*S-\x' l ) (0*eH* n ). 

The vector h n ()3.2|) is nothing but the image of the density matrix p n under the 
latter identification. 

In passing we note that, for a Hopf subalgebra U of U'Js^), A G End(PU) belongs 
to the trivial representation (i.e., &dx(A) = e(x)A for all x G U where e is the 
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counit) if and only if x ■ A — A • x for all x G U. In this case we say A is invariant 
under U. 

With this preparation, let us rewrite our main formula in the matrix formulation. 
Suppose 0* e £* and v G V(Ai) ® • • • V(X n ) ® K(A n + 1) • ■ ■ V(Ai + 1) are identified. 
Then the action v i— > Li(fi)v is translated to the left multiplication 

0* i Li(^)0*, 

while v I— > ^(/i — l)u is translated to 

0* h-> -CTLiOu)- 1 . 

In view of the cyclicity of the trace, the action of the 'transfer matrix' (|3.11j) on v 
turns into 

(-l)"Tr Ml2 T„(/i)(0*), 

where /x = (/ii + /i2)/2 and 

T n (//)(0*) = L n (/2 -\ n )...L x {ii- Ai) ■ 0* • £!(// - AO" 1 • • • L n (ji - An)" 1 . 

Notice that in this formula the normalization factor of the L operator cancels out. 
Regard the operator Q n as acting on £* via the above formula. Then the density 
matrix can be written as 

where I is the identity operator. 

Similarly, denote by the operator corresponding to fl n , acting on the space of 
local operators SL n . Then the main formula (|3.15|) can be rewritten as 

(4.2) (vac|0|vac) = (e n ™(0)) , 
where 

(4.3) n* n (\ u ■ ■ ■ , A n ) = ± J! ^^Tr 2 , 2 (B(^ 2 )(I ® TM)(TM ® /)) 

x (a; 1 (/i 12 )X* in (/ii,/i 2 |Ai, • • • , A n ) + c t ; 2 (/i 12 )X2 )ri (/ii, a* 2 | Ai, • • • , A„)) , 

with 

TM(0) = L 1 (fi-\ 1 )- 1 ---L n (fi-\ n y 1 L n (//-A n )...L 1 (//-A 1 ), 

and X* n , %2 n are constructed from T*(/i) as before. 

Let us discuss the last formula briefly. Consider the operators which act as 
identity either on the last or on the first site (i.e., = 0'£§)/or0 = /£§) 0'). For 
such operators, we have 

(4.4) ST n (I ®0') =10^(0'), 

(4.5) n;(o / ®j) = n^.i(o / )®/- 
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Eq. (|4.5jl is obvious from the definition, while (|4.4jl is non-trivial and follows from 
(jA.lOj) . This motivates us to consider a 'universal' operator 

T»(0) 

= lim L- N (n - A-tv) -1 ■ ■ ■ L N (fi - Ayv) -1 L N (ji - X N ) ■ ■ ■ L- N (p - X- N ). 

N—*oo 

Introducing further the normalized trace 

and defining Q* using T*(/x), we can write down a universal formula 
(4.6) (vac|0|vac) = tr (e n * (0)) , 

which does not refer to the size n of the subsystem. For any operator acting on a 
finite sublattice, the right hand side reduces automatically to this sublattice due to 
(|4.4j) - (|4.5J) . Consider the action of the universal T*(/i) on a given operator 0. If 
acts on sites 1, • • • , n, then the infinite right tail of L-operators Lj with j > n cancels. 
However, the left tail with j < 1 remains. In integrable quantum field theory, 
this situation is typical for the action of non-local charges which transform local 
operators into non-local ones with an infinite tail in one direction [20] • Nevertheless, 
a beautiful feature of our construction is that, when we substitute T*(/i) into the 
trace and integrate, we obtain the operator Q* which sends a local operator to a 
local one because of (|4.5j) . 

In our opinion, this is the most important property of our construction which 
deserves further understanding. 

5. Invariant operators 

As we have seen in the previous section, U'(s\?) operates on our space of local oper- 
ators ()4.ip . The algebra U'Jsh) contains two subalgebras isomorphic to {/^(sfe), one 
generated by eo, fo, q ±h °/ 2 and the other by ei, /i, q ±hl l 2 . In this section, we consider 
the subspace of local operators which are invariant under one of these subalgebras, 
and show that their correlation functions do not contain the transcendental function 
cu 2 (A). To fix the idea, let us choose the subalgebra generated by e , fo, q ±h °l 2 and 
set 

£7 = {0EZ n \x-0 = 0-x (x = e ,fo,q ho/2 )}. 

In the present context, it is more convenient to use the formula for fi* using the 
invariant trace (see the end of Section EJ) , 

(5.1) ^n(Ai, • • • ,X n ) = ^ / / 

x (u;i(/ii 2 )Xi* n (/ii,/! 2 |Ai, • • • , A n ) + W2(/ii2)X2* n (/ii,/i2|Ai, • • • , A n )) . 
As before, X\* n , %l* n are constructed from 7*(/j i ) using Tr^. 
Lemma 5.1. The space £™ v is invariant under the operator fi*(Ai, • • • , A n ). 
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Proof. First we show that Tr*T*(/x) preserves the space £™ v . Abbreviating argu- 
ments, we write 

T n {0) = T- 1 (1®0)T, 

where T — L n - ■ ■ Li, Lj = Lj(fi — Xj) and fi = (/ii + [x 2 )/2. 

Until the end of the proof, we let x stand for eo, fo, q h °^ 2 . The operator T belongs 
to U q (s\ 2 ) ® End(IU) where W = V(Xi) <g> ■ ■ ■ <g> V(A n ). It satisfies the intertwining 
property 

where A(x) = x \ ® x 'l- This equation can be rewritten as 

T(l ®x)= YSrf ® <)T(5- 1 (x:) ® 1), 

where we have set (A ® I) o A(x) = £■ ® (g> a;"'. Using the invariance of and 
the intertwining property again, we obtain, in the notation above, 

(5.2) T-\l ® 0)T(1 ® x) = Y^Wl ® ® ® 1). 

We have q~ h °xq h ° = S 2 (x) and 7T\(q~ h °) = q H , from which follows the invariance 
property 

Tr^(ad'xA) = e(x)Tr^(A) . 
Taking Tr^ of both sides of (|5.2|) and using the above invariance, we find that 

x-Tr*T*(0) =Tr£T*(0) ■ x. 
Let us show that the operator 
(5-3) Tr 2 , 2 (BVi, 2 )(/ ® ^(^ 2 ))(^(Mi) ® /)) 

also preserves £™ v . The relevant operators act on the tensor product C 2 <S> C 2 <S> W . 
To unburden the notation, let us write T 12 = (T(/ii) ® J)(J(8)T(/^ 2 )), = _B 9 (/i 12 ) 
and ©3 = 0, indicating the tensor components by the suffix. Thus the action of 
(E3D on is Tr 2j2 (B^T^O^) . 
Writing again A(x) = J^x^ <S> x", we have 

(5.4) $>^«) 2 ^ = E^(^.)i«) 2 = e(x)Bl 2 . 

Using (|5.4|) together with the intertwining property of T and the invariance of 0, 
we find 

B{ 2 T^0 3 T 12 x 3 = Y,K2^<h T i2^T 12 A{S-\x' l )) l2 . 
Taking trace, moving the last factor by cyclicity and using (J5.4j) again, we obtain 
Tr 2 , 2 (I?Vi, 2 )(/ ® TM){TM ® I)) (0) ■ x 
= ^Tr 2i2 (A(5- 1 (x:)) 12 (xn3^ 2 Tr 2 1 03T 12 ) 
= x ■ Tr 2i2 (^(/i 1)2 )(J (8 TMWM <8> I)) (0), 
which was to be shown. □ 
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Lemma 5.2. For an invariant operator £ £™ v , we /iawe X2* n (j«i, A*2)(0) = 0. 

Proof. Lemma means that the trace 

TrlJlaQi - A^ 1 ■ ■ ■ L n (fi - \ n )- l OL n {^ - A n ) ■ ■ • - A a )) 
= Tr^ i2 (Li(Ai -//)••• L n (A n - /i)OL„(A n - /i)" 1 ■ ■ • Li(X 1 - fi)- 1 ) 

does not produce terms proportional to /ii 2 . We prove this assertion by passing to 
the vector language and performing a gauge transformation. 

Under the isomorphism (|4.1|) . an invariant operator is sent to a vector v £ 
V(Xi) Cg) ■ • • Cg) V(An) ® ^(An — 1) ® • • • Cg) V(Ai — 1) invariant under the action of 



e 



,/ ,g W2 . Set(7 = g 



(1/2) E" =1 (^<+(^-i)^f) 



, and introduce the gauge transformation 



gfog' 



f = ge g 



-i 



q 



h/2 



Q 



-h /2 



£'(X) = q {Xl2)aZ i{X)q 



-(A/2) CT 3 



Then v' = gv belongs to the subspace (y® 2n y uv Q f vectors invariant under U q (sl2 / 



S2nV nv 

lu liic auuopauc i v 

generated by e, /, q h ^ 2 . 

For the proof, we show the following slightly more general statement: for any 
v' £ (V® 2n ) mv and Ai, 



A2n, 



K(A( A l)---4n(A 2 „))^ 



±A2n 



belongs to C[g ±At , g ±Al , ■ • • , q 

First consider the case n — 1, choosing i>' to be s 9 = qv + ®v^ — t>_Cg>t> + £ (V® 2 ) 1 
Note that 



Tr«(F£A) = Tr« 
Tr^EFA) = Tr« 



fJL + l-H' 



l-H 



2 

1 + H 



Direct computation using these relations shows that the entries of ^i(Ai)£ 2 (A2)si 2 
can be reduced to elements of the subalgebra generated by q~ H , Eq~ H ^ 2 , Fq~ H l 2 . 
Since q H does not appear, Tr 9 does not produce terms proportional to \i. 

In the general case, the same argument shows that the assertion holds for v' = 
s q = s\ 2 • • ■ s\ n _ l2n - From the Yang-Baxter relation, the same is true for vectors 
obtained from s q by acting with an arbitrary number of matrices r' ii+1 (Xi i+i), where 

f' (A) = P ii+ iq X ° i/2 r ii+ l{X)q-^' 2 



The operators f' ii+1 (X) are linear combinations of 1 and the generators e, 
(i = 1, • • • , 2n — 1) of the Temperley-Lieb algebra, and vice versa. It is well known 
that the space (y® 2n ) mv is generated from s q by the action of the Temperley-Lieb 
algebra. Hence the assertion is true for all v' £ (y® 2n ) mv . □ 



In summary, let us present the final result for invariant operators. In the notation 
of section HI we have 
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Theorem 5.3. Consider a local operator invariant under the action of U q (s\2) 
generated by eo, fo, q h °^ 2 . For such an operator, we have 

(5.5) (vac\Q\vac) = tr (e n ™(0)) 

where 



(5.6) n; 



1 f f dfii dfj, 



2 



""' 2k 2 J J 2m 2m 

x^(MTr 2 , 2 (B*(^ 2 )7*M ® 7*M) Tr*^ (^1±^1 

Formula (J5.6)) is quite similar to the one in the XXX case jH]. 

Various methods are known for constructing bases of invariant operators (see 
e.g. EI]). For example, a basis for n = 3 is J, Uu, U23, U12U23, U23U12. Here 
the operators Ua + \, after the gauge transformation U' ii+1 = g'Ua + ig _1 by g' = 
qi 1 / 2 ) Ej=i A J CT j ; are the (negative of the) generators of the Temperley-Lieb algebra, 

K +1 = \ + afal, + q -±f^(*W +1 - 1) +^f-(-of + <4i)) • 

(It is nothing but the local density of the Pasquier-Saleur Hamiltonian, see below). 
In the homogeneous case, their expected values are 

(vac|£/j 5 j + i|vac) = — Kao, 

(vac|L r i2?/23|vac) = (vaclL^f^lvac) = (koq — 3k 3 <22), 

COS7T^ 

where k = simru / (my) and 

- tro sinh(l/,-l) trft 
.oo smh(t/z/) cosht 



Let us explain the physical meaning of correlation functions of invariant operators. 
We shall consider the case of massless regime q = e mv (0 < v < 1), because in the 
present context it is more interesting from the point of view of physics. It is well 
known that in the continuous limit the massless XXZ model is described by CFT 
with the central charge c—1. In the continuous field theory, starting with CFT with 
c = 1 one can obtain CFT with c = 1 — by modifying the energy-momentum 
tensor and introducing screening operators. There is a construction which gives a 
lattice version of this procedure, and which is closely related to the invariance under 
the quantum group. 

Consider first the model in the finite volume. The XXZ Hamiltonian with the 
usual periodic boundary condition 



1 L 
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is not invariant under the quantum group. However, following Pasquier and Saleur 
(221 one can introduce an invariant Hamiltonian with specific boundary conditions: 

1 1 

H l V = 2 Yl ( a K+i + + cos(7ri/)<7?<7? +1 ) + -i sm{irv)(al - a 3 ). 

In the infinite volume limit, this Hamiltonian has CFT with c = 1 — as the 
continuous limit. In the finite volume there are significant differences between 
and H™; for example, Bethe Ansatz equations are very different. However, following 
the general logic (see, for example, j2H|) we believe that in the infinite volume the 
ground state of the invariant model is obtained by projection of the original ground 
state onto the invariant subspace: 

|vac)j nv = !Pi nv |vac), 

where Ti nv denotes the projection operator. Then the correlation function of any 
operator in the invariant model coincides with that of an invariant operator in the 
original model: 

inv (vac|0|vac) inv = (vac| (TinvOJmv) |vac). 
So, the correlation functions considered in this section describe the lattice version 
of CFT with c = 1- The fact that they can be expressed in terms of a 

single transcendental function was predicted in [Hj. Certainly, the most interesting 
question is that of rational v when the space of local operators is restricted. We 
shall consider this situation in the future. 

6. XYZ MODEL 

A considerable part of the previous sections can be generalized to the case of the 
XYZ chain 

oo 

(6.1) H XYZ = - (^A+x + + h^l+x) . 

k=— oo 

In jl] we put forward a conjectural formula for the density matrix in the elliptic 
setting. In comparison with the XXX and XXZ chains, however, the results obtained 
are incomplete: we have so far been unable to verify that the expression written down 
in ||] satisfies the reduced qKZ equation. In this section, we rewrite the formula 
conjectured in jl] into a form close to ()3.15|) - ()3.16|) . 

In the following, we denote by 8 a (t) (a = 1, 2, 3, 4) the Jacobi elliptic theta function 
with modulus r. We fix a generic complex number 77 and use the scaled spectral 
parameter t = Xrj. We deal with functions which have period 1 in the variable t. 
The parameters in the Hamiltonian (|6.1|) are given by (1/2)J„ = 9 a +i(2rj) / 9 a +i(0) . 

In the XYZ chain, the role of U q (sl>i) in the XXZ chain is played by the Sklyanin 
algebra (2H1- It is an associative algebra A generated by four symbols S a (a = 
0,1,2,3) satisfying the defining quadratic relations 

(6.2) [5*0, S a ] = Ubc(S b S c + S c S b ), 

(6.3) [S b ,S c ] = i(SoS a + S a S ). 
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Here (a,b,c) runs over cyclic permutations of (1,2,3). The structure constants 
Jbc = —{Jb — Jc)l ' J a are parametrized as 

_ e a+1 (2 V )e a+1 (o) 
a+ M 2 ' 

The algebra A has two basic central elements, K = Yla=o an d K 2 = Yl!=i Ja>3a- 
We will consider only representations of A on which these central elements act as 
scalars, 



(6.4) 

Here and after we set 



2. ' 

[t\- 



d + 1 



d 



The parameter d plays the role of the 'dimension'. For each non-negative integer k, 
the Sklyanin algebra possesses an analog rr^ of the (k + l)-dimensional irreducible 
representation of s^, on which the above relations are valid with d — k + 1. 
The L operator associated with the XYZ chain has the following form. 

p(t,d) 



L(t) :-- 



£(t) 



' [t + {d/2)r,} 
1 3 

2^ e a+ M 



e a+1 {2t+r)) 



S a 



Here p(t, d) is a normalization factor (see 

p(t,d)p(-t,d) = 1, 

P(t,d) p(t-r),d) 



!3l)) which satisfies 



1 



[t + (d/2)rj) [t-n+ (d/2) V ) [(d/2)rj - t] [(d/2) V + t] ' 

We define the monodromy matrix by the same formula ()3.9|) . 

The R matrix R(t) for the XYZ chain is defined by the above formula with 
d = 2 and S a being represented by a a . We need the following three functions 3 
Ui = u>i(t, rj, t) which enter the formula for h n : 

(6.5) Kdlogip = Uidt + U2dr] + u 3 dr, 

where 

0i(2v) 



29' 1 (0) , 
<p(t) :=p(t,2) 



lv-t] 
iv + t] 



1/4 



In the XYZ case, the trace functional is defined as follows. For each element A of 
the Sklyanin algebra, there exists a unique entire function Ti\A with the properties 

i 



3 Note that uj^s here are those in 0] multiplied by fc/4. 
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(i) TrA^L| A _ d = tTir^(A) holds for all positive integers d, 

(ii) If A is a monomial in the S a of homogeneous degree n, Tr\A has the form 



9A,o(t) ( n: odd 



1 ; 



even 



(6.6) TtiA = e 1 (t) n x 

\9A,i{t) - pA,2{t) ( n: 

where gA,o(t), gA,2(f) and gA^if) '■= gA,i(t + T ) —gA,i{t) are elliptic functions 
with periods 1, r. In addition, gA,i(t + 1) = gA,x{t)- 

Let us introduce X a>n (a = 1,2,3) by 

l^sw/lj \ 2 )^ = -^ 1 i«( Sl ' S 2^1) ' ' ' 5 tn) —^2,n(Si, S2 • • ' , £n); 

-^3,n.( s l> S 2 |tl, • • • , t n ) — Xi >n {Sl + T, 52^1, " " " ,t n ) ~ Xl >n (si, S 2 |tl, - • • ,t n ) , 

where X a>n (a = 1, 2, 3) have the following periodicity. 

X a , n (si + 1, s 2 |ti, • • • , t n ) = X a>n (si, s 2 \h, ■ ■ ■ , t n ) (a = 1, 2, 3), 

-^o,n( s l + r ! s 2^1) ' ' ' , tn) — X a>n {Si, S2\ti, - • • , t n ) (a = 2, 3). 

Proposition 6.1. The conjectural formula in [1] for h n can be written as 
_ , , (— l) n /" /" ds\ ds2 ( / w / 1 

(6.7) ilnih,-- ■ ,t n ) = 2 / / ——(^^(^X^!,^^, 

xTr 2 , 2 (T n (si)®T n (s 2 ).S(si i2 )), 



v a=l 




The rewriting procedure is sketched at the end of Appendix IA. 31 

Appendix A. Connection to previous results 

Here we give the details about the derivation of the formula ()3.16|) and the second 
one using (j3.19|) and ()3.20|) . First, in subsection IA.11 we recall the previous result 
in [3], where an algebraic formula for a solution of the reduced qKZ equation is 
constructed. Next we rewrite the formula into the exponential form in subsection 
IA.2l fsee Theorem IA.3|) . Finally we obtain the integral formula (J3.16)) and the second 
one using (j3.19j) and ()3.2()|) in subsection IA.31 The above rewriting procedure is 
applicable also to the elliptic case. We discuss it briefly at the end of subsection 
31 



A.l. Algebraic construction of a solution to the reduced qKZ equation. 

In [3] a solution {h n }^ =0 of the equations ()3.3|) - ()3.5|) is constructed in an algebraic 
way. Let us recall it here. 
First we define the operator 

«x£§(Ai, . . . , A n ) G End(V^ n - 2 \V^ n ) 



X 



3h 
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for 1 < i < j < n as follows. Set 

(A.l) n X n _ 2 (Al, . . ., Xn)(u) :— TT 7- T- — - 

,n,n 

where 

^(A;Ax,...,A n ) := £ I (A-A 1 -1)---A(A^AT-1)---4(A-A n -1) 

x £ n (\-x n )...e t (x^x t )---e 1 (x-x 1 ). 

Then n-^n-2 i s defined by 
nXn^l(\i, • • • , A„) 

:= R^'^(Ai, . . . , A n ) ■ „X„_ 2 (Aj, Xj, Ai, . . . , Aj, . . . , Aj, . . . , A n ), 

where 

^n J H^l) • • • > A ") := ^i,i-l(Ai,i-l) ■ ■ • Ri,l(K,l) 

X #7 ,3-1 K^j ,3 -l) ' ' " ' ' ' ^,1 

X -Rj3Tj(A;_ M ) ■ ■ • ^-(Ai^) 

and 

:= -Pi,i-1 " " " ^2,1 " -Pj,j-1 ' ' ' ^3,2 • ^1=1,1 ' ■ ' -Pl,2 ' Pj=I3 ' ' " ^2,3" 

Note that the operator Kn'^(Ai, . . . , A n ) is a rational function in Q = e mvX] (j = 
1, . . . , n) because p(A)p(— A) = 1. 

From the definition (J3.1U|) of Tr^, the operator n X^l can be written uniquely in 
the following form: 

(A. 2) n X^l(Xi, . . . , A„) = nG^-lid, ■ ■ ■ , (n) — Ajj • nG^'lllCi, . . . , ( n ), 

where n Gn-\{Ci-, ■ ■ • > Cn) an d n^n-lCCi) • • • > Cn) are rational functions in £i, . . . , £ n . 
Take some meromorphic functions cUj(A) (j = 1,2) and consider the operator 

n^i-2( A l5 • • • 5 A n) := w l(Aj,j) • nC^^Cl; • • • j Cn) + ^(Ajj) • n^Vi^Cl; • • • > Cn)- 

For an ordered set of indices K = {ki, . . . , fc m } (1 < k\ < ■ ■ ■ < k m < n), we use the 
abbreviation 

®>K,(k u kj) = m^m- 2 (^fci, • • • , Afe m ). 

Define 

(A.3) h n {X 1 ,...,\ n ) 

[n/2] (_l)m 

:= 2"-^ y^^i.(H.Ji) ' ^K 2 ,(i 2 ,j 2 ) ^K m ,(i m ,j m )( S n-2m)- 

m=0 
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Here 

m 
p=l 

and the second sum in ()A.3|) is taken over all sequences i±, . . . , i m ,ji, . . . ,j m of 
distinct elements of K\ = {1, . . . , n} such that 

h < ■ ■ ■ < i m , h < ji, ■ ■ , im< jm- 

The sets K 2 , . . . , K m are defined by Ki = \ {ii,ji} inductively. 
Now we state the main results of [3]: 

Theorem A.l. The functions {/i n }^L defined by (|A.3J) satisfy the equations (|3.3j) - 
(|3.5|) if u>i(X) andu; 2 (A) are solutions of the following system of difference equations: 

(A.4) wi(A-l)+a; 1 (A)+pi(A) = 0, 

(A.5) cu 2 (A- 1) + oj 2 (X) +wi(A) +Pi(A) +p 2 (A) = 0, 

(A.6) w x (-A) = wi(A), w 2 (-A) = -w 2 (A), 



where 



MX) , 3 1 1 131 



4[A][A-1] 4[A-2][A + 1]' 

P2(A) := 



2[A-l][A-2] 4[A-1][A + 1]' 

Moreover, h n gives the correlation function (|3.2|) m £/ie massive regime when ujj(X) (j ■ 
1,2) are gwen fry ()3.13j) . 

A. 2. Another formula. In this subsection we construct an operator 

n n {\ u ...,\ n )eEnd{V® 2n ) 

and rewrite the formula (jA.3|) by using Q n (see (jA.ll|) ). The procedure of rewriting 
here is similar to that in Section 11 of [3], but slightly different. 

First note that the vector s x . 2 sj j2 in the right hand side of (|A.1|) can be replaced 
by si )2 si j2 because of the following reason. We have the equality 

S 1,2 S 1,2 = s l,2Sl,2 — Si > JS 2j2 . 

From the crossing symmetry 

and the cyclicity of the trace function, we see that 
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Use the quantum determinant formula 



(A.7) Tr d (x£(X)e(-X)) 
Then we find 

As a result we obtain 



Tr d (x) (Vx G E/,(sZ 2 )). 



Tr Al , 2 (tW(^))^ = 0. 



^(i^^)) W2 = Tr Al , 2 (tW(^±^)) Sl , 2Sl , 2 . 

In the following we use the formula for n ^n-2 where s^si^ is replaced by si,2Si,2- 
Introduce the operator n _ifl n e End(V® 2n , V® 2i ~ n - 1] ) defined by 

^ n,n U = {n—l^Ui u )l ,...,n-l ,n-l S n,n, 

and set „_2lT n := n „ 2 n n „i • „-iII„. 

Now we define the operator X^ j) e End(V® 2n ) by 

X^(A 1; ...,A n ) 



(A.8) 



where 



:= -4 n X { Ji(X u . . . , A n ) • n _ 2 ft n • PM R^ j \X u ■ ■ ■ , A n ), 



~p(i,j) ■— p p p p 

1 n ' n— l.ra— 2 J J n,n—l 1 j 



i+l,i ' r n,n-l ' ' ' r j+l,j ' - r n-2,n-l ' ' ' ^1,1+1 ' ^n-l,n ' ' ' ^j,j+V 



and 



^'^(Ai, . . . , A n ) :— Ri, n (\,n) - - " ' ' ' 

x Rj,n\Xj,n) ' ' ' 

X i^j(A nj ) • • • Rj^j(X j+1J ). 



More explicitly, we have 

(A.9) X%' j \\i, ...,X n ) :— 



^ j) (Ai,...,A„) 



■3,Pi 



■ Xi 



x Tr Al , 3 ( e 3 & - 1) (±±±; Ai, . . . , A n ) 



x 



2 7 n v 2 
n^^ J) (Ai,...,A„), 



where 



^n jl (A; Ai, . . . , A n ) 

:= £i(A - Ai - 1) • • -4(A^aT- 1) • • • £j(X^Xi- 1) • • -4(A - A n - 1) 
x 4(A - A n ) • • • ^(A^Aj) • • • £i(\^~Xi) . . . 4(A - Ai). 
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To obtain (lA~9j) from (TOJ) we used 
for it E V mn . 

In the same way as (jA.2|) the operator X n is decomposed into two parts: 

-^n'^C^lj • • • ) ^n) — -^l'n (Clj • • • ) Cn) — Ay " -^'n (Cl j • • • > Cn), 

where (a = 1, 2) are rational functions in £i,...,£ n which take values in 

End(V r ® 2n ). Then take solutions wi(A) and w 2 (A) of the equations (lA~5j) - (lA~6l) . 
and define 

n(«)(A lt . • • , A n ) := £ ^ a (A M )xS')(Ci, . . . , Cn). 

a=l,2 

By definition we set f2« = fin for 1 < i < j < n. Then the operator fln'^ has 
the following properties: 

Proposition A. 2. (1) Suppose that i < j andk<l. If {i, j} n {k,l} 7^ 0, we have 

n^\\ 1 ,...,\ n )- n n i n k i(\ 1 ,...,\ n ) = o. 

If {i,j} H {A;, /} = 0, we /tcroe 

QP(A 1) ...,A n )- n 0i fe i(A 1 ,...,A„) 

= — 4 n fi n l ^ 2 (Ai, . . . , A n ) • n _2^i_4 (Ai, . . . , Ai, . . . , Aj, . . . , A n ) 

x n-i^-n-2 ■ Pn-V ^1-2^1, • • • , Afe, . . . , A/, . . . , A n ), 

where k' and I' are the positions of X k and Xi in (Ai, . . . , Aj, . . . , Xj, . . . , A n ) ; and i' 
and j' are the positions of Aj and Xj in (Ai, . . . , Ajj., . . . , Aj, . . . , A n ). 

The operators fin satisfy the exchange relations: 

Rk,k+l(Ak,k+l)R~k+l ,fc(Afc+l,fc)^n (■ • • > Afc, Afc+l, • • •) 

= . . , X k+1 , X k: . . .)R k , k+1 (X k , k+1 )R^(X k+ i, k ). 

Here Tr k is the transposition (k, k + 1). 
(3) The operators fin are commutative: 

n<M)fi<J.O = f or a u i < j j k<l. 

The proof of Proposition IA.2I is quite similar to that of Lemma 12.1, Lemma 12.2 
and Lemma 12.3 in In the proof of (1), use the recurrence relation 

(A.10) n-lTI n ■ n X n ^ 2 {Xi, . . . , A n ) 

if j = n, 

n-lX^l^Ai, . . . , A n _i) • „_ 3 IIn-2 if j < n. 
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Now we define the operator fl n by 

n n (Ai,...,A n ) := Yl ^ M) (Ai,...,A„). 

l<i<j<n 

Then from Proposition IA.2I and the equality 

^1,2-^2,1-^1,2 (A)-R 2 j( — A)SiJS 2) 2 = •SijS 2i 2, 

we get 

Theorem A. 3. TTie formula (|A.3|) /or /i n can be rewritten as follows: 
(All) /i n (Ai, . . . , An) = 2- n e B »( Al '~' A *>s n . 

A. 3. Derivation of the new formula. Finally we prove that the operator Q„ 
defined in (j3.16j) is equal to fl n . 

From the definition (j3.1H) of X a>n (a = 1, 2), we can easily see that 

res Atl= A J res (U2= A 3 X 0jn (/ii, /i 2 ; Ai, . . . , A n ) = 

and 

res Ml= A i res At2= A J A ajn (/ii, /i 2 ; Ai, . . . , A n ) 
= (-l)°res /il=A Tes /Lt2=A< X 0)n (/ii, /i 2 ; Ai, . . . , A n ). 

The operator 

Tr 2j2 (T n (iii) ® T n (/i 2 ) • S(/ii j2 )) 
is skew-symmetric with respect to \i\ and /x 2 because of the commutation relation 

Ri 2 {ij)B 12 (ij,) = -B 21 (~ij,)R 12 (ij,). 
The functions u a (a = 1, 2) satisfy the parity conditions 

u a (-X) = (-l) a -W(A). 
By using the above properties we have 

fi„(A 1 ,...,A n )= £ Q^' j \X 1: . . . , A n ), 

l<i<j'<n 

where 

(— l) n 

fi^ j) (Ai, . . ., A n ) := — V, ^(AijOres^^res^^Xa.n^i,^; Ai, . . . , A n ) 

a=l,2 

xTr 2j2 (Tn(Ai)®Tn(A i )-5(A i)i )). 

In the following we show the equality Vtn'^ = fin ■ To this end we prove two 
lemmas. 

Let B be a linear operator acting on C 2 <S> C 2 . Then 

Tr 2>2 (T n (Ai) <g>T n (\j) ■ B) 
is the operator acting on V\ ® • • • <S> V n ®V n ® • • • ®V\. 
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Lemma A. 4. We have 

(A.12) iMTnM® 1 ^).- 8 ) 

= ar^\x u . . . , Aj/i,.;)',;)'^^^,. . . . , A n ). 

Proof. We prepare some notation. Denote by 14 eg) V& (V a ~ H ~ C 2 ) the tensor 
product of the two-dimensional spaces on which the trace Tr2 ; 2 is taken. Set 

V Q = V x ® ■ ■ - ® V n ® V n <S> ■ ■ ■ ® Vj. 

Let i? 0ji (A) be the i? matrix acting on the a-th and the i-th component of V^CgV^®^^- 
We set 

T a (A) = R a ,i(X — Ai — 1) • • • -R a ,n(A — A n — l)i2 ai „(A — A n ) • • • -R a ,i(A — Ai). 

We define T b (X) similarly. We denote by the operator B acting on V a eg) V&. 
Similarly, I?,,;,, etc., are defined. 

Denote by X the left hand side of (jA.12|) . Then we have 

X = tT a)b (B aib T a (Xi)T b (Xj)). 

Here tr a fe means taking trace on V a eg) V b . We use the following properties of the R 



matrix. 






(A.13) 


i?i, 2 (0) = Pi, 2 , 




(A.14) 


i2i, 2 (-i) = -2?r, 2 , 




(A.15) 


-^1,2^1,2)^1,3^1,3)^2,3^2,3) = -^2,3^2,3^ 


l-Rl,3(Al,3)-Rl,2(Al,2 


(A.16) 


i?i,2(A)i? 2 ,i(-A) = 1, 




(A.17) 


Ri,s(X)Ri, 2 (-l) = - i? 3 , 2 (-l - X)RiA- 


-1), 


(A.18) 


Pl, 2 (-l)i?l, 3 (A) = - i?l, 2 (-l)i?3,2(-l " 


-A). 



We abbreviate the arguments of the i?-matrices. They can be understood from 
the space indices a, b, 1, . . . , n, 1, . . . , n through the correspondences a <-> Aj, b «-> Aj, 
1 <-> Ai,. . . ,?2 <-> A n , 1 <-> Ai + 1,. . . ,n A n + 1. For example, by we mean 
Rij(Xi j — 1). We extend this convention to T a and T&. We also use the abbreviation 

Ri,[i-l,l] — Ri,i-lRi,i-2 " ' "-R^l, 

etc.. 

With the above convention, we have 
X = tr a fi{B a ,b 

X Ra,[l,i-l]Ra,i( — ^)Ra,li+lj-l]RajRa,[j+l,n]R^ 

x Rb,[i+l,j-l]Rb,j(~ 1) ^6 ) [i+l,n]^6,[n,j'+l]-P6,j ^b,[j-l,i+l]^b,i^b,[i-l,l]}- 

We write the argument —1 in two places in order to emphasize where we can use 
the crossing symmetries (|A.17[) . (|A.18|) . Using the cyclicity of the trace, we bring 
R a ,[i-i,i\Rb,[i-i,i\ to the left of B a b . Then, we move P a ^ to the left while we change 
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all a to i. Finally, we can eliminate tr a P a ^ because this is equal to the identity 
operator. Thus, we get 

X = i^,[i-i,i]tr 6 {i? 6 ,[i-i,i]-B ii6 

X Ri,[l,i-l]Ri,t(~ ^)Ri,[i+lJ-l]RijRi,[j+l,n]Ri,[n,j+l]Ri,jRi,[j-l,i+l] 

x Rb,[l,i-l]Rb~i Rb,[i+l,j-l]Rb,]( — 1) Rb,[j+l,n}Rb,[n,'j+l} Pb,j Rb,\j-l,i+l]Rb,i} 

Using the crossing symmetries ()A.17j) . (|A.18|I . and pushing the underlined terms to 
the right, we obtain 

X = Ri^i^R^zi^-itibiRbXi-i^Bifi 

X Ri,i(~^) R[n,j+l],iRj ,iR[j-l,i+l],iRi,{n,j+l] 

X R[j-l,i+l],jRiJ R[i-l,l]jRb,j(~^-) R[n,j+l],jRb,{n,j+l]Ri,jRj,[j-l,i+l]Rj,[j-l,i+l]Rj,i Pb,j } 

Using the Yang-Baxter equation ()A.15j) and the unitarity relation (|A.16J) . we reduce 
X = -Rj i [j_i i i]i?[j- 1 j]jtr fe {_R b [j_ ljl ]i?j i;) 

X Ri,i(~ l)R[n,j+l] ,iRi,[n,j+l]R[j-l,i+l],]R[j-l,i+l],i 
XR[i-l,l},jRb,j{. — ^-) R[n,j+l],jRb,[n,j+l] Rj,[j-l,i+l} Ri,[j-l,i+l]Pb,j} 

The underlined terms can be brought outside the trace. Therefore, handling P^j in 
the same way as P a ^, we obtain 

X = Ri^i-i^R^^jRijZi^^jRj^j-ij+^Rj^i-i^Bij 

X Ri'i ( — 1 ) R[n, ,lRi,[n,j+l] R\j-l,i+l],i 
X R[i z l,I],]Rj,j(~l) R[n,j+l],jRj,[n,j+l]Ri,[j-l,i+l] 

This completes the proof. □ 
To state the second lemma we introduce some notation. The ^-operator is written 



as 



where 



3 

i(X) = ^f*(A)®o- Q 
«=o 
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For a, P = 0, 1, 2, 3 we set 



i if (a, /3) = (1,2), (2, 3), (3,1), 
-i if (a,/?) = (2,1), (3, 2), (1,3), 
1 otherwise. 



Lemma A. 5. Let .B G End(C 2 (g> C 2 ) be an operator 

3 

B = ^B a \o a ® a?) (B af3 G C). 

Then we have 



(A. 19) Tr Aij ( t n (^±^) ) ftW^Aj 



, A„ ) Bj a 7 . -{P . -. 

1 ILJ t,J j j 



' :J / 



where Y(X) is given by 
(A.20) 



Y(X) := -4 £ - l)n~ - 1). 

a,/3=0 



/ c oo\ 




A 


-l 


-l 


-i\ 






c u 


1 


l 


-l 


l 


l 






C 22 


~ 2 


l 


l 


-l 


l 




£22 


K° 33 / 




V 


l 


l 


-v 







2 y v 2 

i/ere i/ie coefficients C a/3 are determined from the equality 
(A.21) 



and ones obtained by changing the indices (00,11,22,33) in the both hand sides of 
(TOT!) to (01, 10, 23, 32), (02, 13,20,31) and (03,12,21,30). In particular, if B is 
the identity, the left hand side of (jA.19|) becomes zero. 

Proof. We start from the left hand side of (jA.19j) . The i?-matrices l[ iJ ' go through 
the ^-operators t n : 

(A.22) R^Tr^ (/, £j :) ifrX^ij (» h (» bJ)7^PJ s . 

Now it is easy to see that the above operator becomes zero if B is the identity 
because of the cyclicity of the trace, the crossing symmetry 

£ l(A )T- 2 = -£ 2 (-A- 1)0^, 

and the quantum determinant formula (|A.7|) . 

Let us proceed the calculation of (jA.22|) . Without loss of generality, we consider 
the case i = l,j = 2: 

Z : = Tr Al]2 (^ (Y-i) h (¥-0 4 1,2] (A) 4 PP) *i (¥)) ^1,2^,1^ 
where 

= Ai + A 2 



2,2 ' 
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We rewrite the last part of Z by using 

3 

o,/3=0 

Note that 

P\, 2^1 ,1^2,2 = ^'l,2^'l l 2-^l,2- 

Using the crossing symmetry, we rewrite Z as 

3 



a,/3=0 



x Tr Al , 2 (Y-i) of*i (¥) ?L2 • 4 1 ' 21 (A) " 4 (V) °& ^2)^1,2- 

Then ()A.19|) follows from the quantum determinant formula (|A.7[) . the following 
Proposition IA.6[ and the equality 

a 2 £(\y-a 2 = -£(-X-l). 

□ 

Proposition A. 6. Suppose that m = Y2'L=o m a aa where m a is a scalar with respect 
to the 2- dimensional space on which the Pauli matrices act. We have 



a a m = -Eai^m 1 ■ a 2 )a a + 2m, 
mo 

3 ™ f^a\t 



ma a = —e a a a (a 2 m t ■ a 2 ) + 2m a , 



where m 1 = Yla=o m a( ' a ) t an< ^ 

1 if a = 0; 
-1 otherwise. 

The proof of Proposition IA.6I is straightforward. 



Now let us prove that fin = fin ■ If B is the operator £?(//) defined in (|3.17j) . the 
corresponding operator Y(X) determined by ()A.20|) becomes the identity. Therefore, 
from the formula (1A.9D. we obtain 



Tr w>2 (Tn(^r^)) ■ Tr 2 , 2 (T„(A,) ® T n (A,) • 5(A M )) 

= X^(A 1 ,...,A n ). 

This implies the equality fin' -7 = fi« ^. 

The proof of (j3.19j) . (|3.20|) is similar. It is easy to see that for the operator B q {ji) 
we have Y(X) = q H . Therefore, the right hand side of ()3.16|) where X ajn (/ii, /z 2 ) and 
B(ji) are replaced by X^ n (yUi,/i 2 ) and B q (fj,), respectively, is equal to fi n . 
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Finally we give a sketch of the calculation in the elliptic case. The ^-operator is 
given by 

3 

£(X) = ^2w a (X)S a ® o a . 

a=0 

Here S a (a = 0, 1, 2, 3) are the generators of Sklyanin algebra and 

e a+1 (2t + r]) 



w a (X) 



26 a+1 (ri) 



where t — Xrj. If we set 



40i (2t + 277)0! (2t - 2r7) ^ a+ i (2r7) 



cr 09 a , 



the corresponding operator Y(X) becomes 

v(\) i(_0im_ ( Qijt-toiMt-Ti) _ el(t-2 V ) 



QjiM) f 0i(t-v)e 1 (t + v ) e\{t) 

29 1 (2t - 277)0! (2t- 277)01(277) V ° V 

where .Ko and are the Casimir elements ()6.4|) . Therefore we obtain 

TrA(xF(A)) = ^Tr A (x) 

for any element x of the Sklyanin algebra. This gives the formula (|6.7|) . 

Appendix B. Relation with the vertex operator approach 

Correlation functions of the XXZ model in the massive regime have been studied 
in the framework of representation theory ^T], The description in the present 
paper differs from the above literature by a few minor points. In this appendix, we 
compare the two in some detail. 

Consider an inhomogeneous six-vertex model where an inhomogeneity parameter 
Q is attached to each column j. By correlation functions we mean those of local 
operators on a single row. There are two equivalent formulations depending on 
whether one uses row-to-row or column-to-column transfer matrices. In this paper, 
correlation functions are expressed as expectation values with respect to the ground 
state of row-to-row transfer matrices. The latter, and hence the ground state vector, 
depend on the Q while local operators do not. In jT^j, on the other hand, column-to- 
column transfer matrices are employed. Their ground state vectors are independent 
of Q. The inhomogeneity is encoded rather in local operators, expressed in the form 
of an insertion of half-column transfer matrices. 

Another minor difference between ^2] and the present paper is that the R matrices 
and Hamiltonians are not identical. The parameter q = e mv of the present paper 
and the corresponding parameter qjM in are related by 



qjM = -q. 
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With this identification, the R matrix Rjm{C) and the Hamiltonian Hjm in ^2] are 
related to R(X) (|2.4|) and H X xz (ED by the gauge transformation 

(B.l) R(X) = (a 3 <g> 1)^/(0(1 ® A H xxz = KH JM K-\ 

where K = U j:even al 

Correlation functions in the present paper are related to the mean value of the 
two expectation values with respect to the two vectors |vac)(i) considered in [T2j . 
Taking into account the gauge transformation (jB.l|) . we have 

n 

(B.2) JI(-e i )(vac|(S_ e - 1 , 61 ) 1 • • • J„| vac) 

i=l 

1 - 

= 2 n* - ^ ' (0(^1^ ■ (^-61,61)1 • ■ ■ ( E -e n ,eJn ' ^ _1 |vac) (i) 

i=0,l j=l 

= II ^ II ( _ ^)»( VaC l(£-6-i,ei)l---(-S-6 n ,e„)n|vac) (i) . 

1=0,1,7': even j: odd 

The correlation functions 

(i) (vacK^^Ji • • • (-B e „,e„)n|vac) (i) 

can be constructed in terms of the vertex operators arising from representation 
theory of the quantum affine algebra U- q = LL^sk) as follows (recall that qjM = 
■-<!)■ 

Denote by Aj (i = 0, 1) the fundamental weights of U- q . Let V(Ai) be the irre- 
ducible highest weight module with highest weight Aj, and = V ® C[C,C _1 ] the 
evaluation module in the principal picture. The vertex operator (of type I) is an 
intertwiner 

$(C) : V(Ai) — ►y(A 1 _ i )®v c . 

Define the components $ e (C) (e = ±) by 

$ e (C) : V(A0 — * V(A 1 _ i ), $(C) u = J2 (MO u) ® « e . 

Then 

(B.3) (i) (vac|(^_ eii?1 )i ■ ■ ■ (£_ eni? J n |vac) w 

= x" 1 ^ tv v{Ai) (g 2B(i) $ en (g- 1 Cn) • ■ ■ ^(g^CO^CCi) ■ ■ ■ <MU) , 
where Q = D® = -(A + A x ) + |, and 

1 _ (g 2 ;g 4 )oo 
X (gWW 5 (gW)oo' 
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The vertex operators satisfy the following relations: 

(b.4) $ e2 (c 2 )$ ei (G)= £ ^M(Ci/c 2 )^^ e ' 1 (Ci)^ 2 (C2), 

e[,e' 2 =± 

(b.5) e D ^-* e (c)-r DW = 

(B.6) ^£$_ e (C)$e(gC)=id. 

The basic relations (|3.3j) - (j3.4|) are simple consequences of ()B.3|) and (|B.4|) - (jB.6|) . 

Appendix C. The scalar factors 

We collect here formulas for the scalar factors which enter the definition of the 
L-operators. 

XX Z case 

• Massless regime (0 < v < 1) 

, \ , s S 2 (l - I - \) S 2 (2 - % + \) 
(CI) p A, d) = — -) 1 f ^ £ 

where S 2 (A) = 5 2 (A|2, 1/V) stands for the double sine function. 

• Massive regime [y G «IR>o) 

(C.2) p (A,rf) = -C- W ^ ^ W C, Joo 



{q 2 - d C 2 )oo (q 4 - d (- 2 ) 

ii=o 

XFZ case 



where = UT =0 (l - zq A3 ) 



• Disordered regime (rj, t G z'R, — irj > 0) 

(C.3) P (t, d) = -^ ffllf'f^ 

where j(u) = T(u, 4r?, t) and 

00 1 _ g27ri((i + l)r+(fc+l)(7-w) 

r(w, er, r) := ]^[ - — e27ri y T+A;(T+lt ) 

i,fc=0 

denotes the elliptic Gamma function. 

• Ordered regime (r/, t G M, r/ < 0) 

rf) = e ~(d-Vnt xp'(t',d), 
where p'(t', d) is obtained from (jC.3|) by replacing t, rj, r with 

/ t ,77 , 1 

= - r/'=-, r' = — , 

T T T 
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